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Abstract: The dwarf galaxy Segue 1 is one of the most promising targets for the indi- 
rect detection of dark matter. Here we examine what constraints 9 months of Fermi-L AT 
gamma-ray observations of Segue 1 place upon the Constrained Minimal Supersymmet- 
ric Standard Model (CMSSM), with the lightest neutralino as the dark matter particle. 
We use nested sampling to explore the CMSSM parameter space, simultaneously fitting 
other relevant constraints from accelerator bounds, the relic density, electroweak precision 
observables, the anomalous magnetic moment of the muon and S-physics. We include 
spectral and spatial fits to the Fermi observations, a full treatment of the instrumental 
response and its related uncertainty, and detailed background models. We also perform 
an extrapolation to 5 years of observations, assuming no signal is observed from Segue 1 
in that time. Results marginally disfavour models with low neutralino masses and high 
annihilation cross-sections. Virtually all of these models are however already disfavoured 
by existing experimental or relic density constraints. 

Keywords: dwarfs galaxies, dark matter theory, supersymmetry and cosmology, gamma 
ray theory. 
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1. Introduction 

The identity of dark matter is one of the most compelling problems facing modern physics. 
A wealth of viable theoretical candidates have been put forward (see e.g. [1, 2, 3, 4]), with 
the majority based on extensions to the standard model (SM) of particle physics. One of the 
more durable suggestions is that dark matter consists of weakly-interacting massive parti- 
cles (WIMPs), thermally produced in the early universe and therefore naturally present in 
approximately the right cosmological abundance. Models of supersymmetry (SUSY) where 
i?-parity is conserved provide a prototypical WIMP candidate in the lightest neutralino. 
Low-energy SUSY is also highly attractive because it generically solves the SM hierar- 
chy problem whilst simultaneously providing a favourable framework for gauge-coupling 
unification and elect roweak symmetry breaking [5, 6]. 

Because the neutralino is a Majorana particle, its self-annihilation opens a potential 
channel for discovery via the observation of annihilation products like photons, hadrons 
and leptons. Self-annihilation rates are proportional to the square of the particle density, 
so any environment with a high density of dark matter is a good prospective target. In 
practice, expected backgrounds from different targets strongly influence their suitability for 
such indirect detection. Dwarf spheroidal galaxies have recently emerged as leading targets 
for gamma-ray detection of dark matter [7, 8, 9, 10, 11], thanks to their high mass-to-light 
ratios [7, 12, 13, 14] and small expected astrophysical backgrounds. 

Segue 1 is probably the most promising object in this respect [15, 16], due to its extreme 
dark matter domination (M/L ^ 1320), relative proximity (23kpc) and high latitude [14]. 
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As with other dwarf galaxies, constraining the density profile of dark matter in Segue 1 
is difficult; being small and faint, very few stars are available to act as kinematic tracers 
of the gravitational potential. Its spatial superimposition upon the leading arm of the 
Sagittarius stream [17] complicates matters further, as do the partially degenerate impacts 
of dark matter, bulk rotation and magnetic fields upon the stellar velocity dispersion [18]. 
Indeed, the status of Segue 1 as a dwarf galaxy rather than a star cluster, and therefore 
its domination by dark matter, have been called into question [17, 19]. We will assume 
here that it is indeed a galaxy, an assertion strongly supported by further recent (but as 
yet unpublished) spectroscopic data [20]. These new data should also significantly reduce 
the uncertainty associated with the density profile of dark matter within Segue 1. As of 
the time of writing, the best available estimate of this proffie comes from Markov-Chain 
Monte Carlo (MCMC) scans of halo parameters and corresponding solutions to the Jeans 
equation, based on line-of-sight velocities of 24 stars in Segue 1 [15]. 

The Large Area Telescope (LAT; [21]), aboard the Fermi satellite, is a high-energy, 
pair-conversion gamma-ray space telescope. The LAT is designed to operate predominantly 
in survey mode, and has been doing so since August 4, 2008. With its energy range 
(20MeV to over 300 GeV) and high spatial and spectral resolution {AE/E ^ 12%, point- 
spread function < 0.1° at 100 GeV), the LAT is well-suited to gamma-ray searches for dark 
matter annihilation. A major undertaking within the LAT collaboration has been to try 
to discover or place limits upon theories of dark matter using Fermi observations of Milky 
Way dwarf galaxies and satellites, the Galactic centre, the Galactic halo and extragalactic 
sources [22, 23, 24, 25, 26]. The detector design also facilitates direct observation of cosmic- 
ray electrons [27], another possibly relevant channel for dark matter indirect detection. 

By any measure, SUSY is an extensive and highly developed addition to the SM, giving 
rise to a wealth of potential experimental signatures beyond dark matter. Any explanation 
of dark matter as a neutralino must therefore satisfy a host of other phenomenological 
constraints. Even within the minimal supersymmetric extension of the SM (the MSSM), 
the low-energy phenomenology of the theory strongly depends upon the particular parame- 
terisation employed in the soft SUSY-breaking sector, and the specific values of the chosen 
parameters. Given that the most general soft SUSY-breaking Lagrangian in the MSSM 
has over a hundred free parameters, one must choose some reduced parameterisation in 
order to make any progress in fitting experimental data. One approach is to employ a 
low-energy effective Lagrangian for the soft breaking terms, with various parameters set 
to zero or made equal for computational convenience (and in order to avoid experimental 
constraints on e.g. flavour-changing neutral currents). The alternative is to choose a spe- 
cific breaking scheme, such as gravity mediation in minimal supergravity (mSUGRA)[28], 
gauge mediation (GMSB)[29] or anomaly mediation (AMSB)[30], where a small number 
of breaking parameters are defined and unified at some high energy, and the masses and 
couplings are run down to low energy using the renormalisation group equations (RGEs) 
in order to obtain phenomenological predictions. 

In this paper we focus on the Constrained MSSM (CMSSM) as a convenient example 
of one such high-energy parameterisation. This scheme is defined at the gauge coupling 
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unification scale (~ 10^^ GeV) in terms of 4 free continuous parameters and one sign: 

{mo, mi , Aq, tan /3, sgn (1-1) 
Here mg is the universal scalar mass, mi the gaugino mass parameter, Aq the trilinear 

2 

coupling between Higgs bosons, squarks and sleptons, tan/3 the ratio of vacuum expec- 
tation values of up-type and down-type Higgs bosons, and sgn/x the sign of the Higgs 
mixing parameter in the superpotential. We choose fi to be positive throughout this pa- 
per. The magnitude of /i is set by the requirement that SUSY breaking radiatively induces 
electroweak symmetry breaking; in this sense the CMSSM differs slightly from mSUGRA 
(where electroweak symmetry breaking is not strictly part of the definition and tan/3 is 
swapped for the parameter B - see e.g. [5]), but for nearly all intents and purposes the 
two can be considered equivalent. 

The CMSSM possesses a number of distinct regions where the relic density of the 
lightest neutralino matches the observed dark matter abundance (see e.g. [31, 32] and 
references therein). The majority of the CMSSM parameter space results in too high a relic 
density; regions producing the correct amount of dark matter are those where some channel 
of neutralino destruction is especially efficient. Until recently, most analyses focused on the 
so-called bulk region at low mg and mi , where neutralino annihilation proceeds efficiently 

2 

by exchange of light sleptons. This region is now mostly ruled out by collider limits 
on sparticle masses and difficulty in meeting Higgs mass limits when both mg and mi 

2 

are small. The stau coannihilation region occurs at low mg, where the stau is almost 
degenerate in mass with the lightest neutralino. Here the correct relic density is achieved 
via co-annihilations between the two sparticles rather than any increase in the neutralino 
self-annihilation cross-section. A similar situation occurs in the stop coannihilation region, 
which exists at large negative Aq. The stau coannihilation region is still viable, but stop 
coannihilation is disfavoured by low-energy experiments and Higgs constraints. In the focus 
point region at large mo, the lightest neutralino picks up a significant Higgsino component, 
opening new annihilation channels and boosting certain coannihilations. Finally, small 
'funnels' of parameter space exist where neutralino annihilation can be increased by a 
mass resonance with one of the MSSM Higgs particles (i.e. where the Higgs in question has 
roughly twice the mass of the lightest neutralino). The focus point and funnel regions are 
still allowed by present experimental constraints. 

Scanning of MSSM parameter spaces is nowadays a highly developed art. Starting 
from simple grid and random scans [33, 34, 35, 36, 37] within slices of the mSUGRA pa- 
rameter space, efforts expanded to MCMC searches of the full CMSSM/mSUGRA space 
[38], later also including the most important SM uncertainties [39, 32]. As nested sam- 
pling [40, 41] has come to replace the MCMC as the scanning technique of choice, hope 
has risen that MSSM scans might now finally be globally convergent [42, 43]. New results 
with genetic algorithms [44], however, suggest that current scanning techniques may yet 
have some distance to go in this respect. Some authors have begun to focus on higher- 
dimensional low-energy effective MSSM parameterisations [43, 45, 46], which provide for 
a broader range of phenomenological consequences but are almost impossible to scan ef- 
fectively without sophisticated algorithms and substantial supercomputing resources. Ex- 
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plorations of SUSY-breaking schemes beyond mSUGRA have also been carried out lately 
using similar parameter scans [47, 48, 49, 50], as have investigations of next-to-minimal 
SUGRA [51, 52, 53]. SUSY scans are generally either based on the Bayesian posterior 
probability [39, 42, 47, 43], the direct use of the frequentist likelihood [36, 32, 49] (usually 
by a analysis), or a simple 'in-or-out' approach to individual points being permitted by 
experimental data [33, 34, 35, 38]. 

MSSM scans have thus far focused on the constraints provided by particle experiments 
and the dark matter relic density determined from the microwave background, sometimes to 
produce corresponding predictions for astronomical observations (e.g. [15, 54, 55, 56]). To 
our knowledge, none have so far included actual constraints from searches for annihilating 
dark matter; this is no doubt because such constraints have only recently come within a 
reasonable distance of model predictions. 

In this paper, we include the first 9 months of the search for dark matter annihila- 
tion in Segue 1 with Fermi in explicit CMSSM parameter scans. We use spectrally and 
spatially resolved photon counts observed by the LAT to directly assess the likelihood of 
the different regions in the CMSSM parameter space, then combine these with laboratory 
and cosmological data to perform global fits to the model parameters. We also provide 
a predicted impact on the parameter space after 5 years of observations. In Sect. 2 we 
describe our analysis techniques, before presenting results in Sect. 3 and conclusions in 
Sect. 4. 



2. Analysis 

2.1 Gamma-rays from neutralino annihilation in dwarf galaxies 

The expected differential gamma-ray flux per unit solid angle from a source of neutralino 
annihilations is (see e.g. [57]) 

d$ 1 + dN] 



^ f 

Here is the neutralino mass, BF is the boost factor due to any unresolved substructure 
in the source, / labels different annihilation final states, dNj /dE is the differential photon 
yield from any particular final state, Cf is the cross-section for annihilation into that state, 
V is the relative velocity between neutralinos, and the integral runs over the line of sight to 
the source. In the absence of any bound states (i.e. Sommerfeld enhancements), massive 
neutralinos move so slowly that they can effectively be considered to collide at rest, allowing 
afV to be replaced with the velocity-averaged term in the zero-velocity limit, {afv)^. 

Three main channels contribute to the spectrum of neutralino annihilation. Through 
loop processes, annihilation can proceed directly into two photons [58, 59] 

^ = 25{E-m^), (2.2) 
or into a Z boson and a photon [60] 

^""'-^ -S{E-m^ + ^), (2.3) 



dE ^ Am^ 
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giving a monochromatic gamma-ray line. A hard spectrum can also be produced by the 
so-called internal bremsstrahlung (consisting of final-state radiation and virtual internal 
bremsstrahlung) , generated when a photon is emitted from a virtual particle participating 
in the annihilation diagram [61]. Finally, continuum gamma-rays can be produced by 
annihilation into quarks, leptons and heavy gauge bosons (including the Z from the 
line), which subsequently decay via tt^ to softer photons. The cross-sections and resultant 
spectral yields for each of these processes are directly calculable from the SUSY parameters 
which define a point in e.g. the CMSSM parameter space (after appropriate RGE running). 
We use DarkSUSY [62] for this calculation. 

The integral and boost factor in Eq. 2.1 are determined by the dark matter distribution 
in the astrophysical source. We use the Einasto profile [63] 



to describe the average dark matter content of Segue 1, where n is the Einasto index and 
Tg and are the scale radius and density, respectively. This profile is somewhat more 
conservative than the traditional NFW [64], in the sense that it is less steep in the central 
regions, leading to generally better agreement with observations of various dark matter 
halos [63, 65, 66]. It is also slightly more dense at intermediate radii. The adopted form 
of the density profile actually makes little overall difference to the expected flux. This 
is because in general, dwarf galaxies will appear either as point sources or very close to 
pointlike to the LAT, meaning that observations mostly probe the full halo rather than 
just the central cusp. 

We use the best-fit values of the scale radius and scale density found by Martinez et 
al. [15] in their recent fits to stellar kinematic data (rg = 0.07 kpc, ps = 3.8 GeVcm"^). 
Since Martinez et al. found no preference for a particular Einasto index, we adopt the 
central value considered in their scans, n — 3.3. We note that the fits were not only 
influenced by the kinematic data, but also by a theoretical prior imposed by assuming 
the same correlation between and ps as seen in subhalo populations of theoretical A^- 
body simulations of cold dark matter structure formation. Whilst this presents no real 
problem, it is encouraging to see that additional kinematic data [20] largely dominate 
the prior in more recent fits. The same authors performed an extensive investigation of 
the possible substructure boosts in Segue 1, showing that all BF values between and 
^70 are compatible with kinematic data and small-scale structure predictions within the 
CMSSM, with the most likely value depending strongly on the particular model employed 
for the concentration-mass relation. We therefore employ two indicative values for the 
boost factor: a rather pessimistic case, BF — 1, and an optimistic case, BF — 50. It is 
important to note that BF — has a very low probability in the results of Martinez et al. 

2.2 Observations and instrumental considerations 



We considered photon events observed in a 10 square degree, stereographically-projected 
section of the sky centred on Segue 1 (RA,Dec. = 151.763°, 16.074° [14]). We applied cuts 




(2.4) 
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on event zenith angles {9 < 105°), energies (100 MeV < E < 300 GeV) and identifications 
(only 'diffuse class' events - see [21]). All data were processed using the same reconstruc- 
tion algorithms and instrument response functions (IRFs) as the publicly-released first-year 
data. Counts and corresponding exposures were placed into 64 x 64 spatial and 14 log- 
arithmic energy bins. The resultant energy-integrated map of photon counts is shown in 
Fig. 1. 

The Fermi-LAT IRFs consist of the effective area, point-spread function (PSF) and 
energy dispersion. We factored the effective area of the telescope into our calculations of 
the exposure for each bin of observed photon counts, using the the standard analysis tools 
available from HEASARC, specifically ScienceTools 9.11^. We convolved our modelled 
gamma-ray fluxes with the PSF and energy dispersion of the LAT using the publicly- 
available FortranQO library FLATlib [67], which was designed specifically for performing 
this task quickly enough to be useful in MSSM scans. Full Fermi-LAT IRFs are defined 
not only as a function of photon energy, but impact angle with respect to the telescope 
zenith (and even azimuthal angle, though the dependence is weak). FLATlib achieves its fast 
convolution by averaging the IRFs over impact angles, allowing the integral over the PSF to 
be cast as a true convolution and performed by fast spectral methods. The energy integral 
cannot be performed in a similar way, because all three IRFs remain energy-dependent. 
FLATlib performs this integral explicitly, using a fast importance-sampling technique which 
utilises the rough resemblance of the energy dispersion function to a Gaussian. For the sake 
of computational speed, we truncated the PSF at a width of 3.2° in the scans we present 
in this paper. With the present IRF set, this is well beyond the LAT's 95% containment 
resolution at e.g. 100 GeV 0.3°) or even 1 GeV 2°). 

2.3 Likelihoods from Segue 1 

The expected spatial extent of Segue 1 in the gamma-ray sky, if it shines with dark matter 
annihilation, is comparable to the width of the LAT PSF. This puts Segue 1 on the bor- 
derline between a predicted point source and a predicted extended source. For every set of 
CMSSM parameters, we computed model spectra at each pixel in the inner 6x6 square 
shown in red in Fig. 1. We took care to explicitly integrate the density profile over the 
innermost 2x2 pixels as a whole, so as to correctly capture the contribution of the very 
centre of the galaxy (located at their vertex). We compared the predicted spectra with 
the observed ones in each of the 36 pixels to obtain a likehhood based on 504 data points, 
which we then included in the total likelihood for that point in our CMSSM scan. We 
chose only to include the inner 36 pixels in the CMSSM likelihood simply because these 
are the only pixels where there is a predicted signal at any significant level. 

All modelled spectra explicitly included contributions from gamma-ray lines, internal 
bremsstrahlung and continuum radiation. To properly model the observed event counts 
in the region around Segue 1, we also took the Galactic and isotropic diffuse emissions 
into account. We used a preliminary form of the GAL PROP fit to the emission observed 
by Fermi [68] to describe the Galactic diffuse emission. The contribution of the isotropic 
diffuse emission, presumably originating from extragalactic sources, is much weaker and 
depends on the Galactic diffuse model adopted. To describe this, we adopted an isotropic 
^Available at http: //fermi .gsfc .nasa.gov/ssc/data/analysis/sof tware/. 
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power law model with index —2.1, derived from EGRET observations by Sreekumar et 
al. [69]. The models recommended by the LAT team were updated recently, and released 
to the Fermi Science Support Centre. At the ^50° latitude of Segue 1, the differences 
between the old and new models are not important for this analysis. The normalisations 
for both backgrounds were set to the best-fit values obtained in the preliminary 9-month 
LAT dwarf upper-limit analysis [24, 70], based on the full 10 square degree region of interest 
rather than just the inner 36 pixels included in our likelihoods. No sources were detected 
in this region in the first 9 months of LAT operation. 

Because of the very low statistics ob- 
served in LAT photon counts towards Segue 
1, a estimation of the likelihood is inap- 
propriate in this case. We calculated the 
likelihood using a binned Poissonian mea- "I 
sure a 



n 



(2.5) 



'J- 



or, recast in the more familiar minus log- 
likelihood form (analogous to half the x^), 

(2.6) 

Here nj and 9j are the observed and pre- 
dicted number of counts respectively, in 
the jth bin. This prescription clearly ac- 
counts for statistical errors by definition, 
but including systematic errors is less ob- 
vious. To do so, one can marginalise over 
an assumed probability density function 
(PDF) of a systematic error in a semi Bay- 
esian manner, treating it as a nuisance pa- 




220 

Galactic Longitude (deg) 
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Counts 



Figure 1: Photon counts observed by Fermi 
in the region around Segue 1 during the first 9 
months of LAT operation in all-sky survey mode. 
Counts are integrated over all energies between 
100 MeV and 300 GeV. The red cross shows the 
exact location of the centre of Segue 1, and the 
red box shows the region included in our likelihood 
calculations. 



rameter. If we consider a systematic error 

that has the impact of consistently rescaling the observed number of counts as rtj enj 
(i.e. a constant percentage systematic error |1 — e|), and assume a Gaussian form with 
width (jg for the PDF of 6, the marginalised log-likelihood is (see e.g. [71]) 




(e0j)^^e-^^^exp 



2 \ de 



e ^ exp 



Tin 



' 2 I (7, 



de 



d6 . 



(2.7) 



(2i 
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The integral is only analytically soluble for 9j < cr~^, which is not generally true when 
dealing with small statistics; we performed it numerically for each likelihood evaluation. 

We included estimated systematic errors from the LAT effective area (/) and our 
modelled spectra (r) by combining them in quadrature, i.e. (Je{Ej) — f{EjY + t^- Note 
the explicit energy dependence of /; for the present IRF set, f{Ej) ranges from 10% at 
100 MeV, to 5% at 562 MeV, to 20% at lOGeV. We interpolated between these values lin- 
early, and assumed the edge values outside this range. We tuned the importance sampling 
algorithm used by FLATlib using slower, more accurate standard numerical integration 
schemes, choosing a sampling efficiency for our specific problem that would introduce an 
overall systematic theoretical error r of no more than 5% in the normalisation of flux pre- 
dictions. Other systematic errors are no doubt also present in the theoretical predictions, 
but we expect the term from the fast integration to dominate. 



2.3.1 Extrapolation to 5 years of observations 

To make predictions about the impact of 5 years of LAT observations, we explicitly assume 
that no excess events will have been observed after this time. There is no correct way to 
rescale Poissonian counts to longer timescales, so the Poissonian likelihood above cannot be 
used when extrapolating to longer observing times. We instead set the 'observed' number of 
photons equal to the number predicted by the background model, using rescaled 9-month 
exposures. This prescription also avoids the erroneous shifts which confidence intervals 
based on Poissonian statistics can sometimes experience due to a downward statistical 
fluctuation of the background. In this case the observed counts become a continuous 
instead of a discrete variable, so the problem of small statistics disappears. The appropriate 
likelihood measure is then once more the 



.2 



2 _ (^model, j ^observed, j ; _ V ^3 J ^\ 

2 j j model j' ' observed j' 

where $modei,j sind $observed,j ^trc the predicted and observed fluxes, (Jmodeij' ^nd (Joiyservedj 
are their standard deviations, and Sj is the exposure. The exposure is itself the product 
of the effective area and observing time. The standard deviation of the predicted flux can 
be estimated as simply the product of the predicted flux and the percentage systematic 
theoretical uncertainty r (5% in our case - see above), (Jmodeij = '^^modeij- The standard 
deviation in the observed flux can be estimated from the standard deviation of the observed 
counts Grij , and the uncertainty on the exposure , giving 



^observedj " J ^ j ' (^-lO) 

Since the underlying physical process is still Poissonian, the best estimate of is in fact 
cfrij = Furthermore, since the uncertainty in the observing time is negligible, as^ can 

be estimated as simply the percentage systematic error of the effective area f{Ej) times 




-8- 



the actual exposure, as. = f{Ej)£j. We then have 



observed, J 





(2.11) 



~^ ^ observed, jf(-^j) ' 



(2.12) 



giving 




(2.13) 



j 



'mo del, j" 



model, J 



We hasten to point out that constraints based on this extrapolation are probably 
overly conservative, as we assume the same background rejection, systematic errors and 
background model for both the 9-month analysis and the 5-year extrapolation. Our overall 
understanding of the instrument will improve over time, as will our understanding of the 
background as Fermi accumulates better statistics on the Galactic diffuse and extragalactic 
components, leading to correspondingly better constraints on the annihilation cross-section. 
Kinematic constraints upon the dark matter density profile of Segue 1 will also improve in 
time [20, 25], which may impact constraints on CMSSM parameters. 

2.4 CMSSM scans 

We scanned the CMSSM parameter space using a modified version of SuperBayeS 1.35 
[42], employing the MultiNest [41] nested sampling algorithm with 4000 live points. In the 
plots we show, all parameters except those shown on figure axes have been marginalised 
over in some way. In the case of the frequentist profile likelihood, this is simply a matter 
of maximising the likelihood in the other dimensions of the parameter space. In the case 
of the Bayesian posterior, the total posterior (prior times likelihood) is integrated over 
the other dimensions of the space (for a review see e.g. [72]). Because we are somewhat 
more interested in the prior-independent profile likelihood than the marginalised posterior^, 
we prefer linear priors on the CMSSM parameters because they are flat relative to the 
likelihood, causing the sampling algorithm to proceed strictly according to the frequentist 
likelihood function. The effects of alternative priors have already been discussed in detail 
for previous CMSSM scans [42]. 

We used DarkSUSY 5.04 for the relic density and indirect detection computations. This 
allowed us to calculate internal bremsstrahlung spectra, and improved the continuum spec- 
trum and relic density calculations. We also improved the interface between SuperBayeS 
and DarkSUSY, most notably pertaining to the energies at which some particle masses were 
defined. 

Apart from the Fermi data, the experimental data and nuisance parameters which 
we included in scans were identical to those in [42] and [44]. SM nuisance parameters 
were the top and bottom quark masses and strong and electromagnetic coupling constants. 

^We make the point, however, that both should be considered if one wants to gain as complete a picture 
as possible of the preferred regions in an insufficiently-constrained parameter space like the CMSSM. 
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Experimental data were precision electroweak measurements of SM parameters from the 
Large Electron-Positron collider (LEP), the relic density from 5-year Wilkinson Microwave 
Anisotropy Probe (WMAP) fits (n^Mh'^ = 0.1099 zb 0.0062 [73]), LEP constraints on 
sparticle masses, LEP constraints on the Higgs mass, the anomalous magnetic moment of 
the muon (g — 2), the Bg — Bg mass difference, and branching fractions of rare processes 
h 57, Bu yT~ and B^ [x^ [jT . Details can be found in [42]. 

In our chosen configuration, completing the integration over the LAT IRFs for a given 
point in the CMSSM parameter space required a similar order of magnitude in processing 
time as a relic density calculation. Since the relic density computation is the main bot- 
tleneck in MSSM scans, this meant that scans took roughly twice as much total processor 
time to complete as a standard SuperBayeS run. One advantage of FLATlib, however, is 
that it can employ the multi-threaded version of the FFTW library [74] , allowing the IRF 
integration to be performed with a considerably greater degree of parallelisation than the 
present relic density routines. 

3. Results and Discussion 
3.1 Fits to Fermi data only 

In Fig. 2 we show results of scans where the likelihood function only included Fermi data, 
LEP measurements of nuisance parameters and the requirements of physicality (the ab- 
sence of tachyons, that the neutralino is the lightest SUSY particle, and that electroweak 
symmetry breaking is induced by SUSY breaking). Preferred values of the neutralino 
self-annihilation cross-section and mass are shown for scans including 9 months of data, 
scans including the extrapolation to 5 years of data, and a control case without any Fermi 
data. Preferred regions are also given for both the pessimistic and optimistic boost factors 
discussed in Sect. 2.1. 

One apparent feature of Fig. 2 is the lack of viable models with large annihilation cross 
sections for large neutralino masses. This feature is present simply because the annihilation 
cross section goes as vtl^ ^ causing it to fall off at higher masses. 

Given the absence of any observed signal from Segue 1, Fermi data clearly disfavours 
models with the highest cross-sections and lowest masses. This is expected, since higher 
cross-sections and lower masses lead to a larger predicted signal. That constraints are 
best at lower neutralino masses is also consistent with the falling sensitivity of the LAT 
with energy above about 50 GeV, and the reduced source statistics at higher energies. The 
improvement in constraints when moving from the current 9 months of data to the 5-year 
predictions is also roughly what would be expected from a \ft improvement in sensitivity. 
This shows that the two different likelihood estimators we employ give consistent results 
(we also checked this explicitly for 9 months of data, finding very good agreement). 

Predictably, the adopted boost factor plays a large role in determining the extent of 
constraints brought to bear on the CMSSM by Segue 1. In the most pessimistic scenario, 9 
months of LAT observations have no impact on confidence regions, as all disfavoured cross- 
sections are larger than allowed by physicality arguments. In the most optimistic scenario. 
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the data disfavours all models with cross-sections greater than ^ 3 x 10~^^ cm^ s~^. Im- 
provements in constraints when moving from BF = 1 to BF = 50 are consistent with the 
factor of 51/2 improvement in sensitivity expected from Eq. 2.1, as the most pessimistic 
constraints lie above the extent of contours in the upper middle panel of Fig. 2. Extrapo- 
lating to 5 years of observations, all values above 10~^^ cm^ s~^ would be disfavoured, as 
would a region extending down below 10~^^ cm^ s~^ at the lowest masses. Once again, we 
caution the reader that this extrapolation does not take into account systematic improve- 
ments in the background and dark matter profile modelling after 5 years, nor in the LAT 
reconstruction algorithms (see Sect. 2.3.1). 

For comparison, in Fig. 2 we also show the previously-presented, preliminary 95% 
confidence level upper limit from 9 months of LAT observations [70] . This limit was derived 
assuming annihilation proceeds only into bb pairs. Apart from the obvious difference in 
overall strategy (upper limits from an assumed final state versus inclusion in explicit model 
scans), our analysis differs from the upper limit one in a number of ways. The upper limit 
was derived assuming a point source for Segue 1, whereas we perform spatial fits; the 
upper limit is based upon an NFW rather than Einasto density profile, and does not 
include systematic errors nor a treatment of the energy dispersion. 

Nonetheless, the areas disfavoured in our scans are broadly consistent with the 9 month 
upper limit, a positive comment on the reliability of both analyses. Our corresponding 
exclusions do however occur at somewhat higher cross-sections than in the upper limit 
analysis (i.e. our exclusion region is above both the extent of coloured contours and the 
black line in the upper middle panel of Fig. 2). This is to be expected, as our ability 
to exclude models is degraded relative to the upper limit analysis by properly accounting 
for the systematic error in the effective area. Because this error is energy-dependent, our 
exclusions also have a slightly different energy-dependence than the 95% upper limit. 

It should be noted that the degree of substructure apparent in the confidence regions 
of Fig. 2 is unlikely to be physical, and is indeed probably something of an artefact of 
the scanning technique (i.e. 'scanning noise'). In the absence of any constraint on the 
annihilation cross-section from the relic density, the vast majority of points providing a 
good fit to the included data lie at much lower cross-sections. This prompts the scanning 
algorithm to concentrate its efforts there, leaving the region in which we are most interested 
somewhat poorly sampled. From a Bayesian point of view, one would say that when the 
relic density is not included, this region sits well above the most likely annihilation cross- 
sections in the CMSSM, so is not meant to be very well sampled by the nested sampling 
technique. 

Because only a small number of models are disfavoured by including just Segue data in 
the likelihood function, there is little overall impact on the favoured values of mo, mi , Aq 

2 

and tan/3 beyond what is allowed purely on physicality grounds. We will show confidence 
regions from global fits only for these parameters. 

3.2 Global fits 

In Fig. 3 we show the result of including the relic density constraint from the WMAP 5-year 
data, along with all other experimental bounds. The effect is to favour models populating 
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Figure 2: Neutralino self-annihilation cross-sections in the CMSSM, in the zero-velocity limit. 
Left: with no constraining experimental data except measurements of SM nuisance parameters 
and physicality requirements. Middle: constraints provided by 9 months of Fermi data on Segue 
1, under the most pessimistic (top) and optimistic {bottom) assumptions about the substructure 
boost factor. Right: projected constraints after 5 years of Fermi observations. Colours indicate 
68% (yellow) and 95% (red) confidence regions. The preliminary 95% confidence level upper limit 
on the annihilation cross-section from 9 months of Fermi data, assuming 100% WIMP annihilation 
into bb [70], is given for comparison (black curve). 

two distinct regions: a broad strip around the canonical WIMP annihilation cross-section 
at 3 X 10~^^ cm^ s~^, and a low-mass region at smaller cross-sections, corresponding to 
models where stau co-annihilations reduce the relic density to the observed level despite 
the very low self-annihilation rates. The models disfavoured by Fermi observations of 
Segue 1 in Fig. 2 are here already strongly disfavoured by the relic density constraint, so 
the additional data from the LAT appears to have little impact upon the preferred cross- 
sections and masses. A slight reduction in the profile likelihoods of the lowest mass, highest 
cross-section corner of the preferred region appears to be present in the extrapolation to 5 
years of data. 

The best-fit point is however rather different in the 9-month scan as compared to 
scans without Fermi data, or with 5 years of mock observations (where we assumed that 
no excess above background will be seen in 5 years). In the 9-month scan, the best fit 
occurs in the focus point region, at a high annihilation cross-section and a low neutralino 
mass {{(Tv) = 1.8 x 10~^^ cm^ s"-*^, = 95GeV), whereas the best fits in the other cases 
are for stau coannihilation models. This difference appears to be the result of a very 
small statistical excess above the modelled background in the 9-month data. Because the 
corresponding confidence regions are not substantially altered despite the movement of the 
best fit, the excess would appear to be consistent with observational (statistical) noise. 
Given the range of Fermi 's sensitivity, it is thus not at all surprising that the best-fit 
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Figure 3: Annihilation cross-sections in the CMSSM which fit ah experimental constraints, assum- 
ing the neutralino to be the dominant component of dark matter. Favoured regions are as implied 
by existing experimental data only {left\ and with the addition of 9 months of Segue 1 observations 
by Fermi {middle). We also show the extrapolated impact of a non-observation of Segue 1 after 
5 years {right). Upper plots show profile likelihoods (where yellow and red indicate 68% and 95% 
confidence regions respectively), while lower plots show marginalised posterior PDFs (where solid 
blue contours give 68% and 95% credible regions). Solid dots indicate posterior means, whereas 
crosses indicate best-fit points. 



would appear at this location, falling right on the edge of the instrument's sensitivity. This 
point may however be an interesting one to watch as statistics improve. 

It is instructive to note the difference in how the co-annihilation region is represented 
in Fig. 3 by the profile likelihood and the marginalised posterior. Because the range 
of CMSSM parameters spanned by the co-annihilation region is quite narrow (i.e. fine- 
tuned), the total number of points in this region found by the scans is not particularly 
high, leading to a relatively low posterior PDF. This is despite the fact that very good 
fits can be found with a reasonably broad range of neutralino masses and cross-sections 
in this region, as evidenced by its size in the profile likelihood plots. In this sense, the 
Bayesian posterior PDF can be seen to penalise the co-annihilation region to a certain 
degree for being fine-tuned. Whether this is a desirable characteristic or not is of course 
a matter of opinion. It is, however, important to recognise that such information is only 
accessible by comparing the posterior PDF and the profile likelihood; the information in 
their combination is greater than the sum of the parts. 

A natural question to ask might be whether more interesting constraints could be 
obtained from Segue 1 by allowing the neutralino to be a sub-dominant component of 
dark matter. Unfortunately, this generally does not add a lot to the discussion when 
considering constraints from indirect detection with gamma-rays. Even though the relic 
density is essentially inversely proportional to the annihilation cross-section, in mixed dark 
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Figure 4: Annihilation cross-sections in the CMSSM which fit ah experimental constraints, assum- 
ing a 'maximally dense' dark matter halo profile for Segue 1. In this case, the halo scale radius and 
density were chosen ~2(j away from the best-fit values derived from stellar kinematic data. Here 
we again assume the neutralino to be the dominant component of dark matter. Favoured regions 
are as implied by 9 months of Segue 1 observations by Fermi (left), and extrapolations to 5 years 
of data assuming no signal from Segue 1 {right). Shadings and markings are as per Fig. 3. 



matter scenarios the density of neutralinos in Segue 1 becomes directly proportional to the 
relic density. The expected signal is then increased due to the larger annihilation cross- 
sections permitted by sub-dominant relic densities, but reduced by the reduction in signal 
due to the reduced galactic densities. The net result is a reduction in the expected signal, 
since the flux (Eq. 2.1) depends upon the first power of the annihilation cross-section, but 
the square of the density. Thus for a decrease in the relic density such that fi^ fi^/X, 
the flux is modified as $ ^ — ^/X. The result is that the favoured cross-sections 

move to higher values, but the constraints from Segue 1 move even further, providing less 
constraining power than when the neutralino is assumed to be the only component of dark 
matter. This argument of course may not hold for points in the parameter space where 
the relic density is not strictly inversely proportional to the annihilation cross-section, such 
as strong co- annihilation or resonant annihilation scenarios. The former certainly are not 
probed by the Segue 1 observations in any case, since they lie at very low annihilation cross- 
sections. In principle though, highly fine-tuned points in the latter scenario could slightly 
modify the impact of the Segue constraints in subdominant situations. As discussed below 
however, our scans do not uncover a significant number of models where such a mechanism 
occurs. 

In Fig. 4 we investigate whether variations in the dark matter profile of Segue 1, 
within the errors of Martinez et al. [15], might also produce more interesting constraints. 
Here we again take an Einasto profile (Eq. 2.4), but instead use parameters corresponding 
to the most dark-matter-rich profile allowed at ^2a (rg = 10 pc, ps = 70 GeVcm"^). 
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Figure 5: Preferred CMSSM parameter regions including Fermi-L AT observations of Segue 1 and 
all other observables. Shadings and markings are as per Fig. 3. Preferred regions are very similar 
whether one considers the existing 9 months of LAT data or extrapolates to 5 years of observations. 



The corresponding constraints on annihilation cross-sections are indeed stronger than in 
Fig. 3, but are still largely dominated by the relic density. This is not surprising, as even 
though the scale density is a factor of 18 higher in this case, the smaller scale radius means 
that the higher density occurs at a smaller radius. In this sense the two parameters are 
partially degenerate; because Fermi probes essentially the whole dwarf halo (as Segue 1 
should appear almost as a point source), and the total mass of Segue 1 is not substantially 
altered by the change in halo parameters, the corresponding constraints are not massively 
improved. The constraints coming from 9 months of data can be seen to cluster more tightly 
around the best-fit point at low mass and high cross-section, but not to the point where 
significant parts of the rest of the parameter space are excluded. This is consistent with 
our assertion above that any excess can be explained in terms of statistical fluctuations. 

The preferred CMSSM parameter regions including all constraints are shown in Fig. 5. 
Given the marginal impact of Segue 1 observations on scans including the relic density, it is 
not surprising that the regions are very similar to those shown in [42] , even when using the 
extrapolation to 5 years of observations. The stau co-annihilation region is clearly visible 
at low mo and mi, separated from the 'focus point' region at larger mg. Scans indicate 

2 

that both regions are equally well-favoured, though the co-annihilation region tends to 
return the best-fit point in most cases. The 'bulk' region is mostly disfavoured by relic 
density and LEP constraints [32], but persists at low masses in our scans, overlapping the 
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CO- annihilation region in the m^-mi plane. The high-probability region at low tan/? in 

2 

the Ao-tan (3 plane favoured by the co-annihilation region shows up as a much smoother 
peak in our scans than in some previous works [39, 32, 42]. We suspect that this is due to 
our use of the upgraded version of DarkSUSY for the relic density calculation. 

The 'funnel' region, where resonant annihilation can become important at very low mi , 

2 

does not show up in our scans here. This is unsurprising, as the nested sampling algorithm 
is designed to sample according to the total posterior mass, and the linear prior places a very 
small scanning weight upon such fine-tuned regions at low mass. Nested sampling routines 
only find this region when using logarithmic priors on mg and mi [42], though normal 
MCMC scans can find it a little more easily (e.g. [39, 32]). On the other hand, standard 
MCMCs and nested sampling implemented with logarithmic priors sample the focus point 
region less densely, causing them to sometimes miss the highest-likelihood points important 
for a profile likelihood analysis. These difficulties are typical consequences using scanning 
algorithms designed for Bayesian analyses to compute the frequent ist profile likelihood; a 
more promising path for frequentist scans appears to be to use genetic algorithms [44]. 
Using genetic algorithms, it seems possible to find all regions in a prior-independent way, 
but the ability to effectively map their surroundings and produce reliable confidence regions 
lags behind other techniques. 

Some recent MCMC scans [49, 50] have not found large focus-point regions which fit 
all experimental constraints well, leading the authors to claim that the co-annihilation 
region is favoured by present data. In these cases, the reduced likelihood in the focus point 
region relative to the co-annihilation region was almost entirely due to the fact that it is 
virtually impossible to produce a good fit to the muon g — 2 with large values of mo in the 
CMSSM. Using nested sampling with linear priors however, and the physics and likelihood 
routines within SuperBayes, one can find points in the focus point region where this effect 
is essentially offset by a correspondingly better fit to other observables [42]. 

4. Conclusions 

We have incorporated fits to 9 months of Fermi-LAT observations of the dwarf galaxy 
Segue 1 into explicit global CMSSM parameter scans. We included gamma-ray lines, 
internal bremsstrahlung and secondary decay, as well as detailed characterisations of the 
detector response, its uncertainties and the observed background. We have also presented 
scans illustrating the estimated impact of a non-observation of dark matter annihilation in 
Segue 1 after 5 years of LAT operation. 

The LAT data disfavour a small number of physically- viable CMSSM models with low 
neutralino masses and high annihilation cross-sections, but results depend strongly upon 
the assumed substructure boost factor in Segue 1. Such models are already strongly dis- 
favoured by relic density constraints. Extrapolating to 5 years of operation and assuming 
the most optimistic boost factor presently allowed by astronomical data, the absence of any 
annihilation signal from Segue 1 would disfavour all models with cross sections higher than 
10~^^ cm^ s"-*^, as well as a number at low mass with cross-sections as low as 10~^^ cm^ s"-'^. 
Even at this level however, the CMSSM models disfavoured by Fermi would already be 
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essentially excluded by existing data from the microwave background and terrestrial ex- 
periments. 
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